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of the weak response of neutron matter in the region of low momentum transfer. The proposed 
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describe matter at non vanishing temperature. The results show that interactions lead to a sizable 
enhancement of the neutrino mean free path in cold neutron matter. The dependence of the mean 
free path on temperature and neutrino energy is also analyzed. 
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I. INTRODUCTION 

The understanding of neutrino interactions with dense 
nuclear matter is critical to the description of a num- 
ber of astrophysical processes. The results of many-body 
calculations carried out using the formalism of Corre- 
lated Basis Functions (CBF) [IH1] have clearly shown 
that nucleon-nucleon (NN) correlations strongly affect 
the weak response of nuclear matter. Short-range cor- 
relations bring about a sizable quenching of the relevant 
transition matrix elements, with respect to the predic- 
tions of the Fermi gas model, while long-range correla- 
tions lead to the excitation of collective modes that turn 
out to be dominant at low momentum transfer. In Refs. 
OH] short- and long-range correlations have been treated 
in a consistent fashion, using correlated wave functions 
and the cluster expansion technique to derive effective 
interactions from realistic nuclear hamiltonians. The ap- 
proach based on these effective interactions allows for a 
unified treatment of a variety of properties of nuclear 
matter relevant to astrophysical applications, such as, 
for example, the transport coefficients and the neutrino 
emission rates [SHI]. 

An alternative approach to the study of neutrino in- 
teractions in nuclear matter is founded on the theory of 
normal Fermi liquid, in which the dynamics is described 
by a set of Landau parameters [8-113]. This formalism, 
while being strictly applicable only in the regime where 
the excitations of the system can be described in terms of 
quasiparticles, allows for a consistent treatment of collec- 
tive excitations and can be easily extended to non van- 
ishing temperatures. 

In this paper we use the Landau parameters obtained 
from the matrix elements of the effective interaction of 
Ref. [5J to carry out a study of the weak response of pure 
neutron matter in the region of low momentum transfer. 



Section II A is devoted to a short description of the 
many-body approach employed to derive the effective in- 
teraction, while in Section |IIB| wc outline the main el- 
ements of Landau theory, and report the results of our 



calculation of the Landau parameters. In Section |Hl| the 
internal consistency of our approach is tested by compar- 
ing the values of spin susceptibility and compressibility 
computed using Landau parameters to those obtained 
from the matrix elements of the CBF effective interac- 
tion. The density and spin-density structure functions of 



cold neutron matter are discussed in Section IV while in 



Section [V] we analyze the neutrino mean free path and 
report the results of calculations carried out over a broad 
density range at temperatures up to 10, assuming non de- 
generate neutrinos. Finally, in Section VI we summarize 
our findings and state the conclusions. 



II. FORMALISM 

A. CBF effective interaction 

The formalism of nuclear many-body theory provides 
a consistent framework, suitable for treating the non 
perturbative nature of NN interaction. Within this ap- 
proach, nuclear matter is modeled as a collection of A 
point-like particles, the dynamics of which are dictated 
by the hamiltonian 



H = yJPL + y v , 



(1) 



3>i 



where Pi and m denote the momentum of the i-th nu- 
cleon and its mass, respectively, v# is the NN interaction 
potential and the ellipses refer to the presence of inter- 
actions involving three or more nucleons. 

The NN potential reduces to the Yukawa one-pion 
exchange potential at large distances, while its behavior 
at short and intermediate range is determined by a fit of 
deuteron properties and NN scattering phase shifts. 

Performing perturbative calculations in the basis of 
eigenstates of the non interacting system requires the 
replacement of the bare NN potential with a well be- 
haved effective interaction. This is the foundation of the 
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approach developed by Bruckner, Bethe, and Goldstone 
(for a review, see, e.g., Refs. [HI Q2]), in which Vij is 
replaced by the reaction matrix G, obtained by summing 
up the series of ladder diagrams describing NN scattering 
in the nuclear medium. The same philosophy has been 
followed in more recent studies, such as those carried out 
within the self-consistent Green's function approach |13j . 

Alternately, non perturbative effects can be taken into 
account replacing the states of the non interacting sys- 
tem, i.e. Fermi gas states \n FG ) m the case of uniform 
nuclear matter, with a set of correlated states, defined as 
(see, e.g., Ref. [14]) 



F\n FG ) 



(n FG \FiF\n FG y/ 2 



(2) 



The operator F, embodying the correlation structure in- 
duced by the NN interaction, is written in the form 



F = S]Jf k 



(3) 



where S is the symmetrization operator accounting for 
the fact that, as the operator structure of the two-body 
correlation functions reflects the complexity of the 
NN potential, in general [fij,fik] ^ 0. 

Within the CBF approach the new basis defined by 
Eq.Q is employed to perform perturbative calculations 
with the bare NN potential. However, the same formalism 
can be also exploited to obtain an effective interaction, 
suitable for use with the Fermi gas basis [31 [5] . The CBF 
effective interaction, v B s, is defined by the relation 



(0|0) 



= (0 FG \K + v cS \0 FG ) , 



(4) 



where \0 FG ) and |0) denote the Fermi gas and correlated 
ground state, respectively, H is the nuclear hamiltonian 
and K is the kinetic energy operator. 

In Ref. [5], v e ff has been derived starting from a 
truncated version of the state-of-the-art NN potential re- 
ferred to as Argonne vis [H3 HE] . The effects of three- 
and many-nucleon interactions, which are known to play 
a critical role in determining both the spectra of few- 
nucleon systems [T7] and the saturation properties of 
isospin symmetric nuclear matter [18] . have been also 
included, following the approach originally proposed in 
Ref. [TO]. 

The energy per particle of both symmetric nuclear 
matter and pure neutron matter at temperature T = 0, 
obtained from v c s in the Hartree-Fock approximation 



E 



= K 



x [v(0) ~ v(p, 



(5) 



Pi)] . 



turns out to be in good agreement with the results of 
highly advanced many-body approaches [5 . 



In Eq. ^ , Kq is the energy of the non interacting sys- 
tem, n(pjis the Fermi distribution function, 



)(0) = - I d 3 x (AiAj|u e ff|AiAj) , 



and 



p) = - [ d 3 x e ipx J2 (A^kfflAjA 



(6) 



(7) 



where x is the distance between the two interacting nu- 
cleons and p = A/V , V being the normalization volume, 
is the density. The label Xi denotes the spin-isospin quan- 
tum numbers specifying the state of the i-th particle. 



B. Landau parameters 

The theory of normal Fermi liquid developed by Lan- 
dau in the 1950s [5] H] describes a uniform system of 
interacting fermions near T = 0. 

The basic tenet of Landau's approach is that, as inter- 
actions are adiabatically switched on, the states of the 
non interacting system smoothly evolve into interacting 
states. As a consequence, the low energy excitations of 
the liquid, dubbed quasiparticles, retain some of the es- 
sential properties of the non interacting Fermi gas. 

The distribution of quasiparticles with momentum p 
and spin er (from now on, the isospin degree of freedom 
will be dropped, as our discussion will focus on pure neu- 
tron mater) is described by the Fermi distribution 



Jit* 



(8) 



where /i is the chemical potential and j3 = 1/kgT, ks 
being the Boltzmann constant. 

The excitation of the system is measured by the depar- 
ture 8n prT from the ground state (i.e. T = 0) distribution 



5n T 



(9) 



Note that Eq.([8]) is an implicit equation for n pcr , as 
the quasiparticle energy e p(T depends on the distribution 
through 



c pcr 



per 



In the above equation, /, 



p'er' 
crcr'pp' 



crcr'pp' OTlp^ fy' 



(10) 



describes the interaction 



between two quasiparticles with momentum and spin per 
and p'cr' , respectively, while near the Fermi surface e pa . 
can be written as 



e pa- = tl + V F (p-p F ) 

where the Fermi velocity is defined as 



v F 



PF 

m* 



dp 



(11) 



(12) 



P=PF 
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and m* and pp are the quasiparticle effective mass and 
the Fermi momentum, respectively. 

The quantity fa-cr'pp' incorporates both the spin- 
dependence and the non-central nature of the NN inter- 
action in the neutron-neutron channel. As a consequence, 
it can be cast in the form !20 



fa 



fp P > + g PP 'Cr ■ cr' + h pp >Si2{<\) 



where q = p — p' and 



5i 2 (q) 



,(*-q)(*'-q) 



(13) 



(14) 



The role of additional non-central contributions has been 
analyzed in Ref. |21j . However, their inclusion is not 
expected to significantly affect the results discussed in 
this paper. 

Within the approximations underlying Landau's the- 
ory, all momenta are restricted to the Fermi surface, 
where the quasiparticle concept is well defined. Hence, 
one can set |p| = |p'| = pp, and describe the depen- 
dence on cos£ = p ■ p' /pp by expanding f pp >, g pp > and 
hppi in series of Legendre polynomials. For example, the 
expansion of /„„/ reads 



■fpp' =^2fo p t( C0S Q 



(15) 



t=t) 



Up to quadratic terms in Sn pcr , the excitation energy 
of the system is given by the expression 



era" pp' Sflpa-Snp' cr' •> (16) 



P [fin" 



F Fi F 2 G 







Gi G-2 



0.04 
0.08 
0.12 
0.16 
0.20 
0.24 
0.28 
0.32 



-0.569 -0.319 
-0.481 -0.452 
-0.374 -0.544 
-0.263 -0.605 
-0.120 -0.675 
0.004 -0.726 
0.092 -0.734 
0.184 -0.763 



-0.087 
■0.143 
-0.180 
-0.230 
-0.307 
-0.351 
0.390 
-0.412 



0.755 
0.844 
0.881 
0.914 
0.978 
1.006 
1.027 
1.039 



0.169 
0.142 
0.107 
0.067 
0.027 
-0.017 
-0.052 
-0.088 



0.094 
0.138 
0.150 
0.155 
0.150 
0.136 
0.121 
0.102 



TABLE I: Landau parameters Ft and Ge of pure neutron mat- 
ter obtained from the matrix elements of the CBF effective 
interaction. 



P [fin" 



H Hi 



H 2 



0.04 0.060 0.004 -0.018 

0.08 0.072 0.019 -0.016 

0.12 0.056 0.036 0.002 

0.16 0.046 0.060 0.025 

0.20 0.038 0.072 0.041 

0.24 0.028 0.078 0.052 

0.28 0.025 0.084 0.060 

0.32 0.019 0.090 0.069 



TABLE II: Same as in Table[I] but for the Landau parameters 
He. 



III. EQUILIBRIUM PROPERTIES OF PURE 
NEUTRON MATTER 



showing that the quasiparticle interaction can be identi- 
fied with the the second functional derivative of the total 
energy with respect to Sn pcr 



fcrcr'pp' r 



S 2 E 

' per Slip' (j 



(17) 



The above equation can be exploited to obtain fa-cr'pp/ 
from the energy computed in Hartree-Fock approxima- 
tion with the CBF effective interaction [see Eq. ^] . Com- 
bining the resulting expression, 



faa'pp' = [v(0) - i}(p - p')\ 



(18) 



with Eqs. ([6]), ([7| and ( |L3] ), one can then obtain the 
coefficients fi, gi and hi appearing in the expansions of 
fpp' ' 9pp' *^nd h pp ' . 

The Landau parameters F, G and H are dimensionless 
quantities obtained multiplying f pp i , g pp i and h pp > by the 
density of states at the Fermi surface, N — m*pp/n 2 . 

The values of Ft, Gi and Hi obtained from the matrix 
elements of the CBF effective interaction are listed in 
Table [I] and [H] for I = 0,1,2. Their density-dependence 
is displayed in Figs. [l]and[2j 



Before discussing the application of Landau theory to 
the calculation of the dynamic structure functions, in this 
Section we report the results of the calculations of a vari- 
ety of equilibrium properties of neutron matter at T = 0. 

We will focus on isothermal compressibility, x p , effec- 
tive mass, m*, and magnetic susceptibility, x° \ which, 
in the static limit, can be related to the density, energy- 
density and spin-density responses, respectively. As these 
quantities can be obtained both from matrix elements of 
the CBF effective interaction and from the Landau pa- 
rameters listed in Tables |T] and |Hj the analysis discussed 
in this Section provides a valuable consistency test of our 
approach. 

Compressibility and effective mass can be computed 
from the expressions 



,p — 



dP 



1 

V 
1 de 



P d Pj P=PF 



(19) 
(20) 



where the pressure is defined as P = —dE/dV, with E 
given by Eq.([5|, and e(p) is the single particle spectrum, 
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FIG. 1: Density-dependence of the Landau parameters Ft and 
Ge of pure neutron matter obtained from the matrix elements 
of the CBF effective interaction. 
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FIG. 2: Same as in Fig. [I] but for the Landau parameters 
Hp. 



that can be consistently obtained from the CBF effective 
interaction within the Hartree-Fock approximation [3]. 

The corresponding expressions in terms of Landau pa- 
rameters are 

1 N 



m — mil 



P 2 1 + F ' 
Fi 
3 



(21) 
(22) 



The magnetic susceptibility can also be computed us- 
ing the CBF effective interaction, following the procedure 
described in Ref. |22) . 

The energy of matter with spin-up (spin-down) neu- 
tron densitiy p-j- (pi), that can be obtained through a 
straightforward generalization of Eq.([5|, can be cast in 
the form [T] 



E(p,a) = E {p) + E a {p)a 2 



(23) 



with a = (pf - p^) I p. 

In the presence of a uniform magnetic field B, the 
above equation becomes 



where B denotes the magnitude of the external field, the 
direction of which is chosen as spin quantization axis, and 
p is the magnetic moment of a neutron in free space. 
From the definition of total magnetization, M — p{p^ - 
Pi) = X a B, it follows that, at equilibrium, 



X = Po 



2 E 
da? 



P 



(25) 



a=0 



In the case of spherically symmetric interactions, the 
expression of the magnetic susceptibility obtained from 
Landau's theory reads 



1 



Xo = Mo 



1 + Go 



(26) 



When non-central forces are present, the neutron mag- 
netic moment is a tensor of the form 



- , : '- ( PiPj Sij 



3 J ' 



(27) 



where po is the medium modified neutron magnetic mo- 
ment. Neglecting the contributions involving px, the ex- 
pression of the susceptibility in terms of Landau's param- 
eters is E31 [Ml 



~2 

Po 



1 



1 + Go 



1 1 {H -Hi) 



8 1 + Go l + G 2 /5 



(28) 



E B (p, a) = E(p, a) - ap Q B , 



(24) 



The effect of the tensor interaction turns out to be neg- 
ligibly small. Using the values of Ge and Hg reported in 
Tables [I] and pi} and setting p,§ = po, one finds that x a 
given by the above equation differs from Xo of Eq.(26) 
by ~ 0.01%. 

In Fig. [3] the neutron matter compressibility, effec- 
tive mass and magnetic susceptibility obtained from Lan- 
dau theory (dots) and from the CBF effective interaction 
(solid lines) are compared as a function of density. Note 
that x p an d X CT are normalized to their Fermi gas val- 
ues, while the effective mass is given in units of the bare 
neutron mass. 

It appears that the results obtained from the two 
schemes are in close agreement with one another over the 
whole range of densities, extending to twice the nuclear 
matter equilibrium density, po = 0.16 fm - ' 3 . 

The left panel shows that in the low-density region, 
where the effective interaction is predominantly attrac- 
tive, the compressibility is larger than in the non inter- 
acting Fermi gas. Increasing p, the effective interaction 
changes sign and the picture is reversed. 

On the other hand, the right panel shows that the in- 
clusion of interaction effects leads to a suppression of the 
spin susceptibility, with respect to the Fermi gas value, 
at all densities. This behavior suggests that repulsive 
interactions are dominant in spin-density channel. 

In the left panel of Fig. |3J the magnetic susceptibil- 
ity of neutron matter computed within Landau theory 
is compared to the results of Ref. [55], obtained us- 
ing the Auxiliary Field Diffusion Monte Carlo approach 
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FIG. 3: Left panel: Compressibility of neutron matter, normalized to its Fermi gas value, as function of d ensity in units of 
po = 0.16 fm~ 3 . The dots and the solid line correspond to th e re sults obtained from Landau's theory [Eq.(21 1] and the equation 
of state computed using the CBF effective interaction [Eq.(19l], respectively. Right panel: same as in the left panel, but for 
for the spin susceptibility, \° ■ The dashed line shows the susceptibility obtained from the dynamic spin structure function 
including only the incoherent contribution. The inset of the left panel shows the density dependence of the ratio between 
effective and bare neutron mass. 



and nuclear hamiltonians including the truncated v' e and 
v'g [IS] forms of the Argonne Vig potential [T5], supple- 
mented with the Urbana IX three-nucleon potential [26] . 
In the right panel, we compare the density dependence of 
the compressibility resulting from our calculations to the 
results of the variational calculations of Rcf. [TJ], car- 
ried out using the full Argonne vis NN potential and the 
Urbana IX three-nucleon potential. Our results appear 
to be in reasonable agreement with those obtained from 
highly refined theoretical approaches, the differences at 
large density being likely to be ascribable to the differ- 
ent treatment of three-nucleon forces, which are known 
to play a critical role at p > pq. 



IV. RESPONSE FUNCTIONS AT ZERO 
TEMPERATURE 

In the long wavelength limit, corresponding to |q| — > 0, 
where q is the momentum transfer, the dynamic response 
is determined by the physics of the Fermi surface. 

At low |q|, the relevant microscopic degrees of free- 
dom, i.e. the particle-hole excitations, can be treated 
in a ring-like approximation scheme using the formalism 
of Landau theory [TO]. The effect of long range corre- 
lations, which are known to be dominant in this region, 
can be included in a fully consistent fashion, taking into 
account the incoherent and coherent contributions to the 
response on equal footing. 

In this Section, we will outline the application of Lan- 
dau theory to the calculation of the weak dynamic struc- 
ture function, i.e. the correlation function of the hadronic 
weak current, needed to obtain the neutrino mean free 
path in neutron matter to be discussed in Section |V} 

For a generic operator O q , the correlation function can 



be written as the sum of two contributions, according to 



dt 
2t7 



(0|Ot(i) O t( )|0) 



(29) 



Y,(0\Oi\ph)(ph\O^\0)S(uj - uj n0 ) + S M p , 
pi: 



where w„o = E n — Eq, Eq and E n being the energies 
of the states |0) and |n), respectively. The sum in the 
first term includes all one particle-one hole eigenstates of 
the hamiltonian, while the contribution of more complex 
multipair states is included in Smp |27j . 

The linear response function can be written in the same 
form 



X(q,w) = £<0|Ot|pfc)<pfc|O q |0) 

ph 

2u; n o 



+ Xmp , (30) 



with r) = 0+. The link between Eqs.(29) and (30) is 



provided by the fluctuation-dissipation theorem, stating 
that, at T = 0, 



5(q,w) 



-Im x(q,w). 

7T 



(31) 



The theoretical calculation of the first term of Eq.(30) 



is based on the linearized Landau-Boltzmann transport 
equation in the presence of an external probe transferring 
low momentum. Working within this framework, one can 
obtain the explicit expression of the response to a scalar 
probe, x pp (q,uj), involving the two Landau parameters 
F and Fx [TD] 



Nn 



Q(X) 



V 1 + [F + X 2 



Fi 



l+Fi/3 



|n(A) 



(32) 
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FIG. 4: Left panel: comparison between the spin susceptibility computed within Landau's theory and the corresponding 
results obtained in Ref. [55] using the Auxiliary Field Diffusion Monte Carlo approach. Right panel: comparison between the 
compressibility computed within Landau's theory and the corresponding results obtained in Ref. [18] using the variational 
FHNC-SOC approach. 



In the above equation 



A = 



(33) 



is the speed of the perturbation in units of the Fermi 
velocity, while f2(A) is the Lindhardt (or polarization) 
function of the free Fermi gas at T = ,28] . The ex- 
tension to include more than two Landau parameters is 
straightforward [29] . 

It clearly appears that all interaction effects are de- 
scribed by the "potential" in square bracket, that can 
be written in terms of Landau parameters and is there- 
fore constrained by the static properties of matter. Note 
that for cj/|q| < vp the denominator of Eq.(32) can van- 



ish, thus indicating the presence of a resonance between 
the external perturbation and the single quasi particle 
excitation. An undamped mode can exist in the region 
corresponding to A > 1, where the imaginary part of the 
response reduces to a delta function. 



A. Density and spin responses 

Neutrino-neutron scattering via the neutral weak cur- 
rent can be described within Weinberg-Salam theory in 
the low-energy limit. The starting point is the interaction 
lagrangian density 



G 

£/0) = -J= ^0)jzO)> 



(34) 



where G^~4.55xl0 7 fm 2 is the Fermi coupling con- 
stant, while 



and 



3z = ^ n r{l-C Al 5 )i> n , 



(35) 
(36) 



are the neutrino current and the third component of the 
isospin current, with Ca — L25, respectively. Treating 
the neutrons as non relativistic particles, one can approx- 
imate the components of j% according to 



(37) 



(38) 



where a 1 



1,2,3) are Pauli matrices in spin space. 



From Eqs. ( 34 )-(|38[) , it follows that the time component 



of the neutrino current couples to the time component of 
the neutron current to give rise to density fluctuations, 
while the coupling between the space components of the 
neutrino and neutron currents leads to spin-density fluc- 
tuations. 

Let (ko,k) and (k' ,k') be the initial and final four- 
momenta of the neutrino, respectively, while the corre- 
sponding neutron four- momenta will be denoted (po,p) 
and (pg,p'). In the non relativistic limit, the scattering 
rate, computed using Fermi golden rule reads 



W(n,w) 



GjrP 



c 



{(l + cos6»)S(q,o;) 
k^kj+Eki + (1 - cos 6 ) <%] S±(q,uj) 



A 

iC A Siji ( k 



(k e - k l£ ) 



(39) 



where q = k — k', u — k — k' Q , k — k/|k| and 9 is the 
neutrino scattering angle. In the above equation 



S±(q,u>) = -[S ij (q,u)±S ji (q,Cj)} 



(40) 



are the symmetric and antisymmetric parts of the dy- 
namic spin structure function 



Sij(q,u>)= J ^e^iatiq^aji-q^)), 



(41) 
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FIG. 5: Left panel: spin-density structure function of neutron matter, as function of A, denned in Eq. (331, at different densities. 
The calculations have been carried out including the Landau parameters with I = 0, 1 and 2. For comparison, the thick dashed 



line shows the results of Ref. [10], corresponding to p 
density structure function. 
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FIG. 6: Left panel: spin-density structure function of neutron matter at p = po = 0.16 fm , as function of A. The curves have 
been obtained including only the Landau parameters with £ = 0, 1 or taking into account the contribution associated with 
I — 2. For comparison, the thick dashed line shows the results of Ref. [10]. Right panel: same as in the left panel, but for the 
density structure function. 



where (. . .) denote the ground state expectation value 
and (Ji(k,t) is the Fourier transform of the spin-density 
operator ^cr^n- 

The dynamic density structure function can be written 
as in Eq.(|41~]), replacing cTi(q,t) with p(q,t), the Fourier 
transform of the density operator f/'nV'n 



S (q,u)) 



dt 
2^ 



'■(p(q,t)p(~q,Q)) , (42) 



and can be obtained from the response of Eq. ( 32 1 using 
the fluctuation-dissipation theorem, Eq.(31). 



Finally, in an isotropic system = SijS, and Eq.(39) 



reduces to 



{(1 + cos d)S(q,ij) 

+ C*l(3- cos 6)S(k,u)} 



(43) 



The structure function S(k,cu) can be obtained from the 



spin response yf 7 through Eq.(31|. The calculation of 



the spin response within Landau theory in the presence 
of non central interactions is discussed in Ref. 30 . How- 
ever, for our set of Landau parameters, the calculated 
neutrino cross section turns out to be nearly unaffected 
by the inclusion of the Hg's. As a consequence, the 
expression of x CT<T employed in this work is the same as 
Eq. ( 32 ) , but with the Landau parameters Ft replaced by 



the corresponding Gg. 

In Fig. [5] we report the spin (left panel) and density 
dynamic structure functions of neutron matter at differ- 
ent densities as a function of A. It clearly appears that 
in the case of the spin structure function the single quasi 
particle excitations is depleted, and a zero-sound collec- 
tive mode sticks out. As shown by the difference between 
the solid and dashed lines of Fig [3J the collective mode 



provides a large contribution to the spin susceptibility. 
This is a consequence of the repulsive effective interac- 
tion in the spin channel. On the other hand, in the case 
of the density structure function the incoherent contri- 
bution exhausts the compressibility sum rule |31j 



UJ £ 



(44) 



2po, 



and there is no zero-sound mode up to density 
even when Fq change sign. 

The results displayed in Fig. [5] have been obtained 
including the Landau parameters Ge (left panel) and Fe 
(right panel) with i = 0, 1 and 2. However, as shown in 
Fig. [6j the effect of adding the contributions associated 
with G2 and Fi is very small. 

Note that the above discussion can be readily general- 
ized to the case of non vanishing temperature. At T ^ 0, 
the ground state expectation of Eqs. ( 41 (-(42 ) are re- 



placed by the corresponding ensemble averages. 
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FIG. 7: Density dependence of the mean free path of a non 
degenerate neutrino with an energy ko = 1 MeV in neutron 
matter at T = 0. The density and spin-density structure 
functions have been computed using the Landau parameters 
F t and G e of Table |T] with I = 0, 1 (solid line) and I = 0, 1 
and 2 (thick dashed line). For comparison the dashed line 
shows the mean free path in the free neutron gas. 



V. NEUTRINO MEAN FREE PATH 

The neutrino mean free path, L„, can be obtained from 
the transport equation. From the definition of the relax- 
ation time associated with the distribution function of an 
excited state, r, 



drir 



it follows that 
1 



= V 



Of 



d 3 k' 



' per 
T 



{W(q,u)[l-n(k')] 
+ W(-q,-u)n{kf)}.. 



(45) 



(46) 



where n denotes the neutrino distribution function. 
Using the principle of detailed balance, stating that 



W(q,w) 



W(q, 



(47) 



and the definition of the scattering rate of Eq. ( 43 ) , one 
can rewrite Eq. ( 46 1 in the form 



1 



G\ 



d 3 q 



(2tt) 3 



[(1 + cos 0)5(qr, w) 



+ Cl(3- cos 6)S(q,u)] 



(48) 



Note that in the above equation we set T = and as- 
sumed that neutrinos be non degenerate. 

Neutrino-neutron interactions can excite a one 
particle-one hole pair as well as a collective excitation 
associated with the zero-sound mode. Within the frame- 
work of Landau theory, the zero-sound mode propagates 
undamped if A > 1 [see Eq.(33l], i.e. if the phase ve- 
locity of the mode exceeds the Fermi velocity vp- In 
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FIG. 8: Energy dependence of the mean free path of a non 
degenerate neutrino in neutron matter at different tempera- 
tures. 



order for this requirement to be fulfilled in the den- 
sity (spin-density) channel, the Landau parameters Fo,i 
(Go.i) must satisfy the constraint [55] 



F > 0, F Q > 



Fi 



l + fi/3 



(49) 



and the corresponding relations for Go and G±. The re- 
sults listed in Table|7]show that in the density channel the 
zero-sound suffers from strong Landau damping. Hence, 
the spectrum of density fluctuations arises from single 
pair excitations only. On the other hand, this is not the 
case for the spectrum of spin-density fluctuations. 

The above discussion obviously implies that the inverse 
neutrino mean free path is the sum of two contributions, 
that can be obtained from Eq. ( 48 1 singling out the con- 



tribution of the collective mode to the spin-density struc- 
ture function. 

Figure [7] shows the density dependence of the mean 
free path of a non degenerate neutrino with an energy 
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FIG. 9: Density dependence of the mean free path of a non degenerate neutrino with energy ko = 1 MeV at different 
temperatures. 
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1 MeV in neutron matter at T — 0. The results have 



been obtained from Eq.(48), using the density and spin- 



density structure functions computed using the Landau 
parameters Fe and Gt listed in Table[l]with I = 0, 1 (solid 
line) and I = 0, 1 and 2 (thick dashed line). Comparison 
with the mean free path in a free neutron gas, displayed 
by the dashed line, shows that inclusion of interaction 
effects leads to a large enhancement of L v over the whole 
density range. 

At T ^ the mean free path can still be written as 
in Eq. ( 48 1 , using the appropriate expressions of the dy- 



namic structure functions, related to the corresponding 
response functions through [compare to Eq.(31|] 



S(q,w) 
<S(q,w) = - 



x w (q,^) 



7T 1 



Im f'(q,w), 



(50) 



where x pp an d X aa denote the density response and the 
diagonal component of the spin-density response tensor, 
respectively. 

It has to be kept in mind, however, that the description 
of the neutron matter response discussed in this work 
only applies to the regime in which collisions between 
thermally excited quasiparticles can be neglected, de- 
fined by the requirement u S> tJ" 1 , where t^ 1 ss T 2 /Tp, 
Tp being the Fermi temperature, is the thermal collision 
rate. As the zero-mode frequency is ~ T, the collisionless 
regime corresponds to T -C Tp. 

Figure [8] shows the energy dependence of the mean free 
path in neutron matter at density p — 0.16, correspond- 
ing to the Fermi temperatures Tp = 35 MeV. The upper 
and lower curves have been obtained setting the temper- 
ature to T = and 2 MeV, respectively. 

The dependence of the mean free path of a neutrino 
with energy ka — 1 MeV upon both temperature and 
matter density is illustrated in Fig. [9] Note that, as the 
the density range 0.04 < p < 0.32 fm~ 3 corresponds to 
Fermi temperatures 14 < Tp < 55 MeV, the collisionless 
condition T <C Tp is always satisfied. 



functions and the cluster expansion technique. 

Our estimates of the static properties of neutron mat- 
ter turn out to be in fairly good agreement with the re- 
sults of both Monte Carlo [35] and variational jT5] calcu- 
lations. In the case of spin susceptibility, the difference 
between our results and those reported in Ref. [25], the 
size of which can be used to gauge the theoretical uncer- 
tainty, is significantly smaller than the quenching arising 
from interaction effects. 

The calculated dynamic structure functions include 
contributions from both quasiparticle excitations and the 
collective zero-sound mode. However, density fluctua- 
tions are strongly damped, and the ^-function peak as- 
sociated with the zero-sound mode is only visible in the 
spectrum of spin-density fluctuations. 

The dynamic structure functions have been used to 
obtain the neutrino mean free path in neutron matter, 
which plays an important role in determining neutron 
star evolution. The calculations, performed assuming 
non degenerate neutrinos, have been carried out over a 
broad range of densities and extended to finite tempera- 
tures in the region T < 10 MeV, where the collisionless 
approximation underlying our approach can be safely ap- 
plied. 

The results show that interaction effects give rise to a 
large enhancement of the mean free path, corresponding 
to a suppression of the response, providing a measure of 
the neutrino-neutron cross section. 

The effective interaction of Ref. [5] has been recently 
improved with the inclusion of three-body cluster con- 
tributions, which allows for a more realistic treatment 
of three-nucleon interactions, based on microscopic po- 
tential models [33 . A systematic comparison between 
the results obtained from the correlated Hartee-Fock 
and correlated Tamm-Dancoff schemes developed in Refs. 
(31 [31 [33] and those obtained from Landau theory using 
the same dynamics will help to shed light on a number 
of unresolved issues, such as the role of multipair exci- 
tations. Additional information will also come will come 
from the comparison to the sum rules of the weak re- 
sponse computed using the AFDMC approach [34] . 



VI. CONCLUSIONS 



We have studied neutrino interactions in neutron mat- 
ter within the framework of Landau theory of normal 
Fermi liquids. The values of the Landau parameters em- 
ployed in our calculations have been obtained from the 
matrix elements of the effective interactions recently de- 
rived in Ref. [5] using the formalism of correlated basis 
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